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Wo present a measurement of the shape of the Z/^* boson transverse momentum (qr) distribution 
in pp Z/^* e~^e^ + X events at a center-of-mass energy of 1.96 TeV using 0.98 fb~^ of data 
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collected with the DO detector at the Fermilab Tevatron collider. The data are found to be consistent 
with the resummation prediction at low qr, but above the perturbative QCD calculation in the 
region of qr > 30 GeV/c. Using events with qT < 30 GeV/c, we extract the value of g2, one of the 
non-perturbative parameters for the resummation calculation. Data at large boson rapidity y are 
compared with the prediction of resummation and with alternative models that employ a resummed 
form factor with modifications in the small Bjorken x region of the proton wave function. 

PACS numbers: 13.60.Hb, 13.38.Dg, 13.85.Qk, 13.38.Bx 



A complete understanding of weak vector boson pro- 
duction is essential for maximizing the sensitivity to new 
physics at hadron colliders. Studies of the Z/j* boson 
production play a particularly valuable role in that its 
kinematics can be precisely determined through measure- 
ment of its leptonic decays. Throughout this Letter, we 
use the notation "Z boson" to mean "Z/7* boson", un- 
less specified otherwise. 

Z boson production also serves as an ideal test- 
ing ground for predictions of quantum chromodynam- 
ics (QCD), since the boson's transverse momentum, qt, 
can be measured over a wide range of values and can be 
correlated with its rapidity. At large qr (approximately 
greater than 30 GeV/c), the radiation of a single parton 
with large transverse momentum dominates the cross sec- 
tion, and fixed-order perturbative QCD (pQCD) calcu- 
lations [H, 0]) should yield reliable predictions. At lower 
qt, multiple soft gluon emission can not be neglected, 
and the fixed-order perturbation calculation no longer 
gives accurate results. A soft-gluon resummation tech- 
nique developed by Collins, Soper, and Sterman (CSS) 
3 gives reliable predictions in the low-gr region. A pre- 
scription has been proposed for matching the low- 
and high-(7T regions in order to provide a continuous pre- 
diction for all values of qt- The CSS resummation for- 
malism allows the inclusion of contributions from large 
logarithms of the form ln"(gy/Q^) to all orders of per- 
turbation theory in an effective resummed form factor, 
where represents the invariant mass corresponding 
to the four-momentum transfer. The CSS resummation 
can be done either in impact parameter (b) space or in 
transverse momentum (qt) space. In the case of 6-space 
resummation, this form factor can be parameterized with 
the following non-perturbative function first introduced 
by Brock, Landry, Nadolsky and Yuan (BLNY) 



SNpib,Q^) 
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gigs ln(100xiXj 



(1) 

where Xi and the fractions of the incident hadron 

momenta carried by the colliding partons, Qo is a scale 
typical of the onset of non-perturbative effects, and gi, 
g2 and (73 are phenomenological non-perturbative param- 
eters that must be obtained from fits to the data. The 
Z boson qt distribution at the Fermilab Tevatron is by 
far most sensitive to the value of g2 and quite insensitive 
to the value of 53. Thus a measurement of the Z boson 
qt spectrum can be used to test this formalism and to 



determine the value of (72 • 

Recent studies of data from deep inelastic scattering 
(DIS) experiments @,0] indicate that the resummed form 
factor in the above equation may need to be modified for 
processes involving a small- a; parton in the initial state. 
Ref. 3 indicates how such a modification would influ- 
ence the qt distributions of vector and Higgs bosons pro- 
duced in hadronic collisions. A wider qx distribution is 
predicted for Z bosons with large rapidity (called "small- 
X broadening" ) . Z bosons produced at the Tevatron in 
the rapidity range 2 < < 3 probe processes involving 
a parton with 0.002 < a; < 0.006, and can be used to test 
the modified form factor at small x. 

Z boson qx distributions have been published previ- 
ously by the CDF [9] and DO collaborations using 
about 100 pb~^ of data at ^/s — 1.8 TeV. In this Letter, 
we report a new measurement with larger statistics and 
improved precision. This measurement is also the first 
to present a qx distribution for large-rapidity Z bosons. 
The data sample used in this measurement was collected 
using a set of inclusive single-electron triggers with the 
DO detector at the Fermilab Tevatron collider, and 
the integrated luminosity is 980 ± 60 pb~^ [l^ . 

Our selection criteria for Z bosons require two iso- 
lated electromagnetic clusters that have a shower shape 
consistent with that of an electron. Electron candidates 
are required to have transverse momentum greater than 
25 GeV/c. The electron pairs must have a reconstructed 
invariant mass 70 < M(ee) < 110 GeV/c^. If an event 
has both its candidate electrons in the central calorime- 
ter (CC events), each electron must be spatially matched 
to a reconstructed track. Because the tracking efficiency 
decreases with rapidity in the endcap region, events with 
one or two endcap calorimeter electron candidates (CE 
and EE events, respectively) are required to have at least 
one electron with a matching track. After these require- 
ments, 23,959 CC, 30,344 CE, and 9,598 EE events are 
selected; 5412 of these have a reconstructed Z boson with 

|y| >2. 

Electron identification efficiencies are measured using 
a combination of data and a GEANT-based [31 simulation 
of the DO detector. The electron identification efficien- 
cies are measured from Z data. The dependence of the 
overall selection efficiency on the Z boson qx is param- 
eterized from the GEANT simulation. A measurement of 
this shape from the data agrees well with the simulation 
within statistical uncertainties. 
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The dominant backgrounds are from photon plus jet 
events and di-jet events, with photons and jets misiden- 
tified as electrons. The kinematic properties of these 
events are obtained from events that satisfy most of the 
Z selection criteria, but fail the electron shower shape 
requirement. The normalization of the background is 
obtained by fitting to a sum of a signal shape obtained 
from a parameterized simulation of the detector response 
and the invariant mass distribution from the background 
sample to the invariant mass distribution of the data 
sample. The background fractions are (1.30±0.14)%, 
(8.55±0.26)%, and (4.71±0.30)% for CC, CE, and EE 
events respectively. Other backgrounds are negligible. 

The data are corrected for acceptances within a range 
of generated Z masses of 40 to 200 GeV/c^, and for se- 
lection efficiencies using a parameterized simulation. We 
use ResBos 1J| as the event generator which incorpo- 
rates the resummation calculation in 6-space using the 
BLNY parameterization for low and a NLO pQCD 
calculation for high qr- We use PHOTOS [Hj to simulate 
the effects of final state photon radiation. The overall 
acceptance times efficiency falls slowly from a value of 
0.27 at low qT to a minimum of 0.19 at qT = 40 GeV/c 
and slowly increases for larger qT- 

The measured spectrum is further corrected for de- 
tector resolution effects using the RUN (Regularized Un- 
folding) program [3| to obtain the true differential cross 
section. Its performance was verified by comparing the 
true and unfolded spectrum generated using pseudo- 
experiments. The measured Z qx resolution is about 2 
GeV/c; the bin width we choose is 2.5 GeV/c for qT < 30 
GeV/c. The typical correlation between adjacent bins is 
around 30%. Due to limited statistics, the chosen bin 
width is 10 GeV/c for 30 < < 100 GeV/c and 40 
GeV/c for 100 < qr < 260 GeV/c. 

Systematic uncertainties on the unfolded qT spectrum 
arise from uncertainties on the electron energy calibra- 
tion, the electron energy resolution, the dependence of 
the overall selection efficiency on qT, and the effect of 
parton distribution functions (PDFs) on the acceptance. 
The uncertainties on the unfolded spectrum are esti- 
mated from the resulting change when the smearing pa- 
rameters are varied within their uncertainties. CTEQ 
6.1M is used as the default PDF. Uncertainties due to 
the PDFs are estimated using the procedure described in 
Ref. 17|. The uncertainty due to the choice of unfold- 
ing parameters in the RUN program is also estimated and 
included in the final systematic uncertainty. 

The final results in the qT < 30 GeV/c range, are 
shown in Fig. [1] for the inclusive sample and for the sam- 
ple with \y\ > 2. Each data point is plotted at the 
average value of the expected distribution over the bin 
[lit . For the theoretical calculation, we use ResBos with 
published values of the non-perturbative parameters [1]. 
Good agreement between data and the prediction is ob- 
served for all rapidity ranges, which indicates that the 



BLNY parameterization works well for the low qx region. 

Z boson events produced at large rapidities {\y\ > 2) 
are also used to test the small- a; prediction. We compare 
data with the theoretical predictions with and without 
the form factor as modified from studies of small-x DIS 
data Q. All curves are normalized to 1 for qT < 30 
GeV/c. The default values for the parameters gi, g2, 
and gs obtained from large-x data are used. The 
X^/d.o.f. between the data and the ResBos calculation 
using the default parameters is 0.8/1 for qT < 5 GeV/c 
and 11.1/11 for qT < 30 GeV/c, while that for the modi- 
fied calculation is 5.7/1 for qT < 5 GeV/c and 31.9/11 for 
qT < 30 GeV/c. It remains to be seen if retuning of the 
non-perturbative parameters could improve the agree- 
ment for the modified calculations. 

Figure [2] shows the measured differential cross sec- 
tion in the range qT < 260 GeV/c compared to (1) the 
ResBos calculation with its default parameters [5|, (2) 
ResBos with a NLO to NNLO K-factor by Arnold and 
Reno tl9|] incorporated into ResBos by its authors, (3) 
a pQCD calculation at NNLO using the MRST 2001 
NNLO PDF set divided by the NNLO calculation of 
the total cross section [2l| , and (4) the NNLO calculation 
but rescaled to the data at qT = 30 GeV/c. The agree- 
ment between data and ResBos, with or without the 
K-factor, is good for qT < 30 GeV/c. At higher values 
of qT, the data are not in agreement with the ResBos 
calculation. The data agree better with the NNLO cal- 
culation and ResBos prediction with the Arnold-Reno 
K-factor, but agrees best when the NNLO results are 
rescaled by a factor of 1.25 so that they match the data at 
qT — 30 GeV/c. This indicates that the shape from these 
calculations agrees with the data, and that the source of 
the discrepancy is in the normalization. Table U summa- 
rizes the measured values for each qT bin together with 
statistical and systematic uncertainties. 

The CSS model parameter most sensitive to the shape 
at low qT {qT < 30 GeV/c) is g2- In a fit, we fix other 
phenomenological parameters to the values obtained in 
Ref. [H| and only vary g2. A minimum x^/d.o.f. of 
9/11 between the model and the inclusive data for qT < 
30 GeV/c is found when 52 = 0.77 ± 0.06 (GeV/c)^. 

In conclusion, we have measured the normalized dif- 
ferential spectrum, for Z boson events produced 
in pp collisions at ^/s = 1.96 TeV with boson mass 
40 < M < 200 GeV/c2 and qT < 260 GeV/c. This rep- 
resents the highest center-of-mass energy measurement 
of this quantity over the largest phase space available to 
date. The overall uncertainty of this measurement has 
been reduced compared with the previous measurements. 
We find that for qT < 30 GeV/c, the CSS resummation 
model used in ResBos describes the data very well at all 
rapidities. Our data with |?/| > 2 disfavor a variant of this 
model that incorporates an additional small-x form factor 
when 51, (72, and g^ from large- a; data is used. Using the 
BLNY parameterization for events with qT < 30 GeV/c, 
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FIG. 1: The normalized difTerential cross section as a function of qr for (a) the inclusive sample and (b) the sample with 
Z boson \y\ > 2 with qr < 30 GeV/c. The points are the data, the solid curve is the ResBos prediction, and the dashed line 
in (b) is the prediction from the form factor modified after studies of small-x DIS data. 
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FIG. 2: The normalized differential cross section as a function of qr compared to four theoretical calculations for (a) the 
entire range measured and (b) the fractional differences between data and the theoretical predictions. The four theoretical 
calculations are ResBos with its default parameters, ResBos with a NLO to NNLO K-factor by Arnold and Reno, the NNLO 
calculation by Melnikov and Petriello, and the NNLO calculation but rescaled to data at qr ~ 30 GeV/c. 



we obtain g2 = 0.77 ± 0.06 (GeV/c)^, which is compara- 
ble with the current world average value ^ . We observe 
a disagreement between our data and NNLO calculations 
in the region qT > 30 GeV/c, where our distribution is 
higher than predicted by a factor of 1.25. However, the 
NNLO calculation agrees in shape with our data when 
normalized at qx — 30 GeV/c. 
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TABLE I: The normalized differential cross section for 
Z events produced in bins of qr- The first uncertainty is 
statistical and the second is systematic. 
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